Triplet-singlet conversion in ultracold Cs 2 and production of ground state molecules 
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We propose a process to convert ultracold metastable Cs2 molecules in their lowest triplet state 
into (singlet) ground state molecules in their lowest vibrational levels. Molecules are first pumped 
into an excited triplet state, and the triplet-singlet conversion is facilitated by a two-step spontaneous 
decay through the coupled A X E+ ~ b 3 H u states. Using spectroscopic data and accurate quantum 
chemistry calculations for Cs2 potential curves and transition dipole moments, we show that this 
process has a high rate and competes favorably with the single-photon decay back to the lowest 
triplet state. In addition, we demonstrate that this conversion process represents a loss channel 
for vibrational cooling of metastable triplet molecules, preventing an efficient optical pumping cycle 
down to low vibrational levels. 

PACS numbers: 32.80.Pj, 33.20.-t, 34.20.-b 



I. INTRODUCTION 

Research on cold and ultracold molecules currently at- 
tracts considerable interest due to their features and ap- 
plications in fundamental science [l| , from the control of 
low-energy quantum dynamics of few-body systems, to 
the intrinsic many-body nature of quantum degenerate 
gases [2J , as well as due to theprospects in applied areas 
such as quantum computing 3] , precision measurements 
[1, Q , and ultracold controlled chemistry [(| . Therefore 
work on creating samples of cold molecules has rapidly 
progressed, and recent reviews outline the experimental 
methods, theoretical issues and app lications for neutral 
P, 0, [H and charged species Except in remark- 

able cases [HI, [l3[ , the rich internal structure of molecules 
prevent them to be directly cooled down by lasers. 

Among the main challenges is the formation of ultra- 
cold molecules in a well-defined internal quantum state, 
which would open the way for the full control of its evolu- 
tion under various interactions. Therefore work on creat- 
ing samples of ultracold molecules has rapidly progressed 
in two main directions: either molecules are created by 
association of ultracold atom pairs induced by a pho- 
ton (photoassociation, or PA) JUfl or by a magnetic field 
(magnetoassociation, or MA) [llj, [l6| , or by direct cool- 
ing applied to pre-existing molecules, such as buffer gas 
cooling [17[ or Stark deceleration [18j. The former ap- 
proach yields molecular samples with temperatures in 
the micro- or nanokelvin range, but most often with high 
internal vibrational energy. In contrast the latter ap- 
proach produces molecules in their lowest energy levels, 
but their translational motion is characterized by a tem- 
perature larger than 1 millikelvin in most cases. The 
ultimate goal to produce ultracold molecules in micro- 
or nano-Kelvin range with no vibrational or rotational 
excitation, or in a well-defined internal quantum state 
still poses a difficult challenge. Using stimulated Raman 
adiabatic passage (STIRAP) technique, the population 
of high-lying bound levels of molecules created by MA of 
ultracold atoms in a degenerate or near-degenerate quan- 



tum gas has been transferred down to the lowest molec- 
ular bound level through a single STIRAP step in KRb 
19] and through a double STIRAP step in Cs 2 @,[H1. 
These results pave the path towards the achievement of 
a degenerate gas of molecules with no internal energy. 
Additionally, ultracold alkali-metal diatomic molecules 
produced by PA have been transferred into their lowest 
vibrational level v = by stimulated emission pump- 
ing for RbCs [lU and LiCs [23| and by optical pumping 
for Cs2 [24j . However the relatively small fraction (a few 
thousands) of ground state molecules in the ultracold gas 
could be a limitation for further studies. 

In all the cases above, coupling between excited elec- 
tronic state mediating the population transfer plays a 
crucial role, as it has been previously demonstrated in 
several PA experiments [25l - l28l |. and some other possi- 
bilities are still to be discovered [29] . Typically, samples 
of ultracold molecules are initially created in their low- 
est metastable triplet electronic state, and are converted 
into ground state molecules via a transfer mechanism re- 
lying on a triplet-singlet coupling in excited electronic 
states induced by spin-orbit interaction. In the present 
paper, we propose a triplet-singlet conversion mechanism 
in Cs2 which takes advantage of the large formation rate 
of triplet molecules that are indeed produced by PA of 
ultracold cesium atoms (30l . l3l| . A one-photon transition 
excites molecules in the lowest a 3 £+ state into a well- 
defined level of the excited 2 3 n s (6s + 5c?) state which 
relaxes down into low vibrational levels of the ground 
state in two steps via the 0+(A 1 E+ ~ b 3 U u ) coupled 
states (FigHJ). We show that the probability for accumu- 
lating molecules in their ground X 1 !]^ state is compa- 
rable to the one for accumulating of molecules back into 
the a 3 E+ state. Besides the triplet-singlet conversion of 
cold molecules, this mechanism allows for u — g symmetry 
conversion, which cannot be achieved via a two-photon 
process in homonuclear molecules, in contrast with het- 
cronuclcar molecules. As a further result, this conversion 
process suggests a likely explanation for the suppression 
of vibrational cooling of a 3 £+ molecules down to low vi- 
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FIG. 1. (Color on line) Scheme of the proposed process, (a) 
the a 3 £j v = 17 level is initially populated by PA. (c) It can 
be efficiently excited (thick vertical arrow) to the (2) 3 n 9 v = 
12 level, (b) the ,4^+ (6s + 6p) and 6 3 n„(6s + 6p) potential 
curves, and the resulting 0^ curves (dot-dashed lines) after 
diagonalization of the potential matrix of the A and b curves 
coupled by the spin-orbit coupling function of ref . [33i] . The 
(2) 3 n 9 level can decay either directly to the a state (arrow 
I), or to the X state via the 0„ states (arrow II and III). 
Probability densities are reported with an arbitrary scale. For 
the coupled wave function, the component on the b state 
(upper red trace) and on the A state (lower black trace) of 
the wave function are both represented on the same scale, (d) 
independent dipole moment for the transitions a 3 E^-(2) 3 II fl 
(D (2) . a ), 6 3 n„-(2) 3 n 9 (D (2) _ b ), and X^j-A 1 ^ (D A . X ). 



brational levels as reported in ref. [321] . 

The proposed model is described in Section Inland the 
transfer efficiency is evaluated in Section HIIl The compe- 
tition with other possible decay channels is discussed in 
Section IIVI where prospects for experimental realization 
of the proposed scheme are considered. 



II. THE PROPOSED TRIPLET-SINGLET 
CONVERSION SCHEME 

We investigate the efficiency of the process depicted 
in FigfTJ Our model is based on accurate spectroscopic 



information recently obtained on Cs2 electronic states. 
The X 1 E+ ground state potential curve is taken from 
ref. while the lowest a 3 E+ potential curve comes 
from the recent analysis of Li Li's group [35j . The triplet- 
singlet conversion is mediated by the A x S+(6s + 6p) and 
b 3 U u (6s + 6p) states (hereafter referred to as the A and 
b states, respectively) coupled by a spin-orbit (SO) in- 
teraction, which depends on the internuclear distance R 
[36| . According to Hund's case c labeling, these coupled 
states result in a pair of 0+ states dissociating into the 
6 2 S'i/2 + 6 2 Pi/2,3/2 limits. In the following we will refer 
to the 0^(A ~ b) pair of states (see dotted lines in FiglTJ . 
The potentials curves for the A and b states are calcu- 
lated by ab initio methods described in ref. [13], which 
were adjusted empirically to reproduce the spectroscopic 
data of references [H[ and [35[. The SO couplin g fun ction 
is taken from the ab initio determination of ref. [331 ] . The 
(2) 3 n g potential curve as well as the independent transi- 
tion dipole moments are also evaluated according to the 
method of ref. [13 (Fig{TJi). The (2) 3 II g potential curve is 
found with an harmonic constant w e = 17.1 cm -1 in good 
agreement with the measurement or ref.[39j; it is shifted 
downwards by 112 cm -1 to match its minimum with the 
value reported in ref. 39] for the (2) 3 ILj(l g ) state. The vi- 
brational energies and wave functions for all these molec- 
ular states are computed with the Mapped Fourier Grid 
Representation (MFGR) method [iof . 

We assume that Cs2 molecules are initially created by 
PA of cold Cs atoms into the so-called giant Gl or G3 
resonances (at 11 720 and 11 715 cm -1 ) characterized 
in ref-HHIll of the double-well 0"(6 2 S' + 6 2 P 3/2 ) state. 
These photoassociated molecules are stabilized by spon- 
taneous emission in a distribution of vibrational levels of 
the a 3 £+(6s + 6s) (hereafter referred to as the a state) 
located around the if —17 level, as discussed in ref. p2l] . 
These molecules are efficiently excited towards low-lying 
vibrational levels of the 2 3 II g (6s + 5c?) state. Starting 
from one of the most populated (v a = 17) level of the a 
state, the strongest transition is found towards the v = 12 
level of the (2) 3 LT ff state at 710 nm or 14084.6 cm" 1 . This 
is illustrated in Fig. [2^ where the squared matrix element 



= 17|L> ( 2)_ a |(2) 3 n g ,'y)| 2 is shown. Then two 



electronic relaxations are open: either a direct transi- 
tion back to the a state, or a two-step transition down 
to the X state via the A — b coupled states. The ar- 
rows illustrating these transitions in FigfT] suggest that 
the classical turning points of the potentials match well 
enough to ensure that these two ways compete together, 
as confirmed in the next sections. 



III. THE TWO-STEP SPONTANEOUS DECAY 
OF THE 2 3 n 9 STATE 

In this section, we demonstrate that the v = 12 level 
of the 2 3 II S state can efficiently create ground state 
molecules in their lowest vibrational level if = 0. We 
consider the problem only from the point of view of the 
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transitions between vibrational levels, where transition 
strengths depend only on vibrational wave functions. We 
left apart the contribution of the rotational state of the 
molecules, which could be incorporated through Honl- 
London factors, without significantly modifying the main 
conclusions of the paper. The decay rate of a vibrational 
level Vi of an electronic state A towards all the vibra- 
tional levels Vj of an electronic state A' is given by the 
general expression for the Einstein coefficient A v . (A — A') 
(in s- 1 ): 



-4, 



,(a-a')=x; 



16tt 3 3 



(i) 



where hwy is the energy difference between the V{ and Vj 
levels, and Da-A' (R) the i?-dependent transition dipole 
moment between the electronic states A and A'. 

We first compute the relaxation rate of a given (2) 3 il ff 
vibrational level down to the 0^(A ~ b) levels. This rate 
obviously depends on the amount of b electronic charac- 
ter for each of the bound levels of the + (A ~ b) cou- 
pled states, and on the dipole moment function D b _( 2 ) 
(FigHH) for the transition (2) 3 n g -6 3 n u . We see in Fig. 
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FIG. 3. Relaxation rate A v (in s _1 ) of the 2 3 Il ff vibrational 
levels (a) towards the levels of 0£(A — b) system, and (b) 
towards the levels of the a 3 T,t state. 



displayed in FigQJi. The envelope of this function reflects 
the oscillatory structure of the v = 12 wave function. The 
v' = 71 level of the coupled system corresponds to a local 
maximum in a range which is predicted below to decay 
favorably down to the v" = level of the X state. This 
is also illustrated in Fig[T] where it is shown that one of 
the turning points of the v' = 71 wavefunction indeed co- 
incides with the minimum of the X potential curve. We 
note that the double STIRAP scheme of ref. [21| actually 
used the 0+(?/ = 63)-X(v" = 0) for their final step of the 
transfer. 

Using an expression similar to EqO] we display in 
Fig. [5^ the decay rate of the 0+ levels down to the 
ground state. It depends on the squared matrix element 
|(A 1 S+w"| J D j4 _x|0+,w'}| 2 - Below v' = 42, i.e. below the 
crossing between the A and b potential curves, 0+ levels 
have mostly a b 3 H u character, so that their decay rate to- 
wards X levels is very low. Above this region vibrational 
levels acquire a significant A 1 !!^ character whose ampli- 
tude oscillates from one level to another (see for instance 
ref. [3(| ) . Such levels contribute to the dark background 
of lines in Fig. [5^,. Among them, the levels which have 
the largest overlap with X wave functions (like v" = 71) 
are visible as lines with an amplitude about two times 
larger than the one of the background, and 4 to 5 times 



larger than that of the (2) II S 



12-0+ (A - b) decay 



Hi that (2) 3 n s 



12 level has a relaxation rate of rate. Figure [Sb shows that the squared matrix element 



9 x 10 6 s _1 over the entire 0+ system. The population 
of the 0+ levels v' after this relaxation is related to the 
squared matrix elements |(0 + , v'\D^ 2 )-b\(2) 3 Ilg, v = 12)| 2 



71)| is one of the largest 



|(X 1 E+ «» =0\D A _ x \0t 
one. 

To summarize, this study shows that we can find an 




FIG. 4. Squared matrix element ^j 2 = 

\{0t , i/|D(2)_f,|(2) 3 n g , v) 2 (in atomic units) of the transition 
dipole moment function Dt-(2) (see Fig[TJl). The level 
v' — 71 discussed in the text is indicated for clarity, (a) 
v = 12; (b) v = 21. 



efficient transfer of the triplet molecules initially created 
by PA down to ground state molecules in their lowest 
vibrational level v" = 0, relying on their excitation into 
a specific level of the (2) 3 Hg state. Our model involves 
the v' = 12 level, but this could slightly vary depending 
on the availability of new spectroscopic investigations, 
without changing the main conclusion of the proposal. 
Then the two-step spontaneous decay down to the v n = 
levels looks promisingly efficient, despite its apparent 
complexity induced by the strong SO coupling in the Cs2 
molecule. 



IV. DISCUSSION 

As shown in Fig[TJ the single-step decay of the (2) 3 II 9 
levels back to a 3 £+ levels obviously competes with the 
above conversion process. The rate depending on the 
squared matrix element (a 3 S+, v a |Z?(2)-al(2) 3 n s , w)| 2 is 
evaluated according to Eq. [TJ and is represented in 
Fig|3p. This rate is only twice as large as the (2) 3 II s -0+ 
rate, which confirms that both processes will indeed take 
place. The v = 12 level (Fig(B^) then decays significantly 
back to the v a = 17 level, which could probably be partly 
rcpumped during the conversion process, depending on 
the details of the experimental procedure. Also, we have 
not considered in our model the \ u and 2 U components 
of the b fine structure manifold for the spontaneous de- 



FIG. 5. (a) Relaxation rates A v > (in s _1 ) of the vibrational 
levels of the 0t coupled states towards the levels of the X 
ground state. The v' = 71 level is indicated for clarity, (b) 
The squared matrix element {(X^+v" = 0\D A -x\0t ,v')\ 2 , 
indicating that the v' = 71 Oj level indeed decay efficiently 
to v" = 0. 

cay of the (2) 3 n s state used for the conversion. These 
states represent a statistical weight of 80% within this 
manifold, which may reduce the computed conversion ef- 
ficiency. On the other hand, the only possible channel 
for these states to decay is the X state, either through 
their admixture of (l) 1 ]!^ state (for the l u component), 
or through weaker couplings (for the 2„ state). 

This conversion process is in principle valid even if we 
start from high-lying vibrational levels of the a state, 
obtained for instance when the PA laser is not exciting 
the specific giant resonances of the 0~(6 2 S , +6 2 P 3 / 2 ) state 
of Cs2, or when molecules are formed by MA. This is 
illustrated in FigJ^b where a quite low level (v = 21) of 
the (2) 3 II S can be populated from the highest vibrational 
v a = 54 level in the a state. However the transition 
probability (proportional to the squared matrix element 
represented in FigJ5]) is about 100 times smaller than for 
reaching the v = 12 level from low lying a levels, due to 
the large radial extension of the v a = 54 wave function 
with a weak amplitude in the short range. The other 
steps of the process are comparable to the previous case. 
Figure 0Jd shows the distributions of the populated 0+ 
levels, which again reflects the oscillatory structure of the 
upper v — 21 function. The v — 71 level is about twice 
less populated than from the v = 12 level. As expected 
finally, the one-step decay of the v = 21 level (Fig. [5]) 
populates individual a levels with a probability about 5 
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transition dipole moment function D a _( 2 ) (see FigUJi). (a) 
v = 12; (b) v = 21. 



light (as in ref. |43|) can cover v a = 15 — 25 levels with 
the most of the molecule population. This process would 
transfer more molecules into the 0+ system effectively in- 
creasing the triplet-singlet conversion and accumulation 
of molecules in the ground X state. 

Additionally, the proposed triplet-singlet conversion 
probably explains the suppression of vibrational cooling 
of the a state [32]. As the (2) 3 II g state was used for 
the cooling transition, our work shows that a large frac- 
tion of the excited molecules are most likely lost at every 
absorption-emission cycle of the cooling process due to 
the two-step decay to the X state. 

The proposed conversion scheme is currently imple- 
mented experimentally at Laboratoire Aime Cotton. The 
transitions involved in this conversion are within an easy 
reach for laser light. Moreover, the outlined process re- 
lies on spontaneous emission, but larger rates could be 
obtained with coherent transfer such as STIRAP, while 
molecules may be prepared in a well-defined internal 
ground state. This kind of conversion process would ac- 
tually create a molecular cold gas suitable for the study 
of cold collisions between atoms and molecules, or among 
molecules, prepared in a well-defined internal state. As 
the number of cold molecules is quite large after the ini- 
tial PA step [3l|, accumulating them in the v = level 
of the ground state would also provide an ideal starting 
point for the study of cold atom- molecules photoassocia- 
tion into cesium trimcrs. 



to 10 times smaller than the v — 12 decay , but generates 
a broader distribution of a levels. 

To minimize the effect of relaxation back to the ground 
triplet a state, molecules can be optically pumped into 
the 0+ system by applying additional laser excitation to 
the a 3 S+ state. A setup employing a broadband laser 
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